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Abstract. The single-site two-electron exchange amplitude Jgd between the 
Cu 4s and Cu 3d^2_y2 states is found to be the pairing mechanism of high-Tc 
overdoped cuprates. The noninteracting paxt of the Hamiltonian spans the copper 
Cu 4s, Cu 3d^2_y2 and oxygen O 2px and O 2py states. Within the standard 
BCS treatment an explicit expression for the momentum dependence of the gap 
Ap is derived and shown to fit the angle-resolved photoemission spectroscopy 
(ARPES) data. The basic thermodynamic and electrodynamic properties of 
the model [specific heat C(T), London penetration depth A(T)] are analytically 
derived. These are directly applicable to cuprates without complicating structural 
accessories (chains, double Cu02 planes, etc.). We advocate that the pairing 
mechanism of overdoped and underdoped cuprates is the same, as Tc displays 
smooth doping dependence. Thus, a long-standing puzzle in physics is possibly 
solved. 



PACS numbers: 74.20.Z, 74.20.Fg, 74.72.-h 

1. Introduction 

The discovery of high-temperature superconductivity [1] in cuprates and the 
subsequent "research rush" have led to the appearance of about 100 000 papers to 
date [2]. VirtuaUy every fundamental process known in condensed matter physics 
was probed as a possible mechanism of this phenomenon. Nevertheless, none of 
the theoretical efforts resulted in a coherent picture [2]. For the conventional 
superconductors the mechanism was known to be the interaction between electrons 
and crystal-lattice vibrations, but the development of its theory lagged behind the 
experimental findings. The case of cuprate high-Tc superconductivity appears to be 
the opposite: we do not convincingly know which mechanism is to be incorporated in 
the traditional Bardeen-Cooper-Schrieffer (BCS) theory [3]. Thus the path to high-Tc 
superconductivity in cuprates. perhaps carefully hidden or well-forgotten, has turned 
into one of the long-standing mysteries in physical science. 

Features of the electronic spectrum of the Cu02 plane, figure 1(a), the structural 
detail responsible for the superconductivity of the cuprates, have become accessible 
from the angle-resolved photoemission spectroscopy (ARPES) [4,5]. Thus, any theory 
which pretends to explain the cuprate superconductivity is bound to include these 
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features and account for them consistently. A number of extensive reviews over the 
past years have been devoted to that theoretical problem [6-17]. For further related 
discussion we also refer the reader to the review [18] on NMR-NQR spectroscopies in 
high-Tc superconductors. 

In contrast with all previous proposals, we have advanced in reference [19] the 
intra-atomic exchange [20] of two electrons between the 4s and 3da;2_j^2 states of the Cu 
atom as the origin of high-Tc superconductivity in the layered cuprates and have shown 
that the basic spectroscopic and thermodynamic experiments can be explained by it. 
Previously only inter-atomic Heitler-London-type [21] two-electron exchange [22-24] 
has been discussed. Thus, the present work is the unabridged version of our theory 
announced in reference [19]. It builds upon the standard Bloch-Hiickel [25-28] (tight- 
binding) approximation to the electronic band structure of the Cu02 plane, developed 
in an earlier work [29]. We derive an analytical expression for the BCS kernel, 
or pairing potential Vppi . For the case of the s-d pairing the analytical solution is 
compared to the ARPES data. Extensive discussion is also provided to help the 
juxtaposition of our theory with other models. Finally, exact expressions within the 
s-d model are given for the specific heat, London penetration depth. Cooper-pair 
effective mass and Hall constant of the vortex-free Meissner-Ochsenfeld phase. 

2. Lattice Hamiltonian 

The electronic properties of materials are strongly influenced by the local environment 
and in this sense the electronic features are local physics. The simplest possible model 
for high-Tc superconductivity contains single-particle and correlated two-electron 
hoppings between nearest neighbours and next-nearest neighbours. Formally, this 
is an expansion of the many-particle Hamiltonian containing two- and four-fermion 
operators. The two-fermion Hamiltonian determines the band structure, briefly 
considered in subsection 2.1, while the four-fermion terms (subsection 2.2) determine 
the pairing interaction, and lead to the gap equations considered in section 3. 

2.1. The four-band model in a nutshell 

Every high-Tc superconductor has its specific properties. It is strongly believed, 
however, that the main features of the electronic band structure of the Cu02 plane 
are adequately described by the four-band model spanning the Cu 3dx2-y2, Cu 4s, 
O 2px and O 2py orbitals, figure 1(b). In the spirit of the Bloch-Hiickel (BH) 
model, using Jordan's second quantization language, we introduce Fermi annihilation 
operators for an electron with spin projection a at a particular orbital, respectively, 
Dna, Sna, Xna, and Yna in tlic Unit ccU with index n = {nx,ny). It is convenient 
to introduce a multicomponent Fermi creation operator in momentum space, = 
{D'^pa, S^pa^ ^pa^^pa) ■ ^his notation the one-electron BH Hamiltonian reads 




(2.1) 



where is the chemical potential, and (cf. reference [29]) 





tpdSy \ 



(2.2) 
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Figure 1. (a) Ball-and-stick model of the Cu02 plane. The shaded square is 
the unit cell indexed by n = {nx,ny), nx,y = 0, ±1, ±2, . . .. (b) The LCAO basis 
set: A single electron hops from the 3d atomic orbital to 2px with amplitude 
tpd, contained in Wbh- From 2px to 2py the hopping amplitude is tpp, and from 
there to 4s the hopping amplitude is tup ■ Correlated hopping of two electrons in 
opposite directions between 3d and 4s with amplitude J^j is depicted as a double 
arrow (see the discussion in sections 4 and 5). 



Cd, Cs, and ep are the single-site energies of the Cu ZA,j;2_y2^ Cu 4s, O 2p^ and O 2pj, 
states, respectively. The hopping amplitudes between these states are tg-p, tpd and tpp, 
figure 1(b). Note, that because of the orbital orthogonality tsd = 0. For brevity, we 
have introduced also the notation 

Sa; = 2sin(p3;/2), Sy = 2Sl\i{Py/2), S^{Sx,Sy), 

X =sin2(p,/2), y^&in\py/2). (2.3) 

From a classical point of view, the Cu 'iA^2_y2 state corresponds to a circular 
electron rotation in the Cu02 plane while the Cu 4s state corresponds to a classical 
ensemble of electrons of zero angular momentum continuously falling to the nucleus. 
Pictorially, the s-electrons fall to the nuclei like comets, but after the impact the 
turning point of their motion is very far from the nucleus. This is the reason why tgp 
is considerably larger than tpd. The transfer amplitude ipp is the smallest one since the 
hopping to the next-nearest neighbour requires a tunnelling through free space. As a 
rule, the electron band calculations significantly overestimate tpp, but the latter can be 
reliably calculated using the surface integral method, cf. reference [29]. Even for the 
largest transfer integrals tgp and <pd, which determine the bandwidth of the conduction 
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band, the ab initio calculations give a factor 2 or even 3 "overbinding". Nonetheless, 
the band calculations substantiate this choice for the LCAO (linear combination of 
atomic orbitals) basis set and provide an adequate language for interpretation. In the 
end, these parameters should be determined by fitting to the spectroscopy data and 
be treated in the model lattice Hamiltonian as phenomenological parameters of the 
microscopic many-body theory. We shall briefly recall some basic properties of the 
four-band model as derived in reference [29]. 

Let ef,,p and '^b,p be the eigenvalues and the corresponding eigenvectors of the 
BH Hamiltonian, i?BH^6.p = £b.p^b,p, where b = 1,...,4 is the band index. For 
the "standard model", ep < ed < Es, the lowest energy band, 6 = 1, is an oxygen 
bonding band having a minimum at the (tt, tt) point. The next band, 6 = 2, is 
a narrow "nonbonding" oxygen band with an exactly (within the framework of the 
model) zero dispersion along the (0, 0)-(7r,0) direction, i.e., this band is characterized 
by an extended Van Hove singularity. The conduction band, 6 3, is a nearly half- 
filled Cu 3dj^2_j^2 band with the analytical eigenvector 







5*3, p 




X3.P 









-Es el + 4£p t^p (x + y)- 32<pp r^^ xy \ 



y) 

-(£s Ep - 8T^py)tpd 

(Es Ep - x) tpd Sy 



(2.4) 



where the e's denote the energies measured relative to their respective atomic levels: 
Es = e — Es, £p = e — gp, £d = e — fid, and t^^ — t'^^ ~ estpp/2. The topmost band, 
6 = 4, is an empty Cu 4s band. In elemental metals like Cu and Fe, the 4s band is a 
wide conduction band, but for the Cu02 plane it is completely "oxidized". Having the 
analytical eigenvector we can calculate the corresponding eigenvalue: 



(^3,p|gBH|^3,p) 



(2.5) 



'3,p|*3,p) 

If necessary, the nonorthogonality of the atomic orbitals at neighbouring atoms can 
be easily taken into account. In this case the normalizing denominator in the above 
equation reads (for arbitrary band index) 

(*p |*p) = Dl + Sl + X^ + Yl + 2gpAS,DpXp - ^g^ASyDpY^ 

"^dsp^xSpXp -\- "^(jsp^y SpYp "ZgppSxSyJCpYp^ (2.6) 

where the "metric tensor" gij is given by the integral 



3u 



^*{r)il;j{r — R)dr, 



(2.7) 



where '>p*{r) and 4'j{r — R) are the atomic wave functions, and R is the inter-atomic 
distance. The phases are chosen such that all overlap integrals ppd, gsp, and gpp 
be positive parameters, like the hopping integrals ipd, tgp, and tpp. Note that these 
provisions apply only to the single-particle spectrum. As long as one deals with a 
single conduction band, all Bloch states are orthogonal and the further treatment of 
the second-quantized Hamiltonian proceeds in the standard way. 

Thus, using the Rayleigh quotient iteration for equations (2.2)-(2.5) one can 
obtain numerically the eigenvalue and the eigenvector. The band energies e = eb.p 
satisfy the secular equation 



det {Hbh - el4x4) = Axy + B{x + y) + C = 0, 



(2.8) 
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where the energy-dependent coefBcicnts read [29]: 

^(e) = 16 (4tpdtsp + 2igpipped - 2tpjjtpp£s - tppedSs) , 



(2.9) 



C(6) 



Furthermore wc mtroducc the normahzed eigenvector '^b,p = ^'fc.p/||5'fc.p|| and write 
the noninteracting Hamihonian in diagonal form, 



BH 



b.p, a 



(2.10) 



The Fermi operators in real-space representation can be easily expressed using the 
band representation. 



\y I 

\ ^ no. / 



7Np 

b.p 



( Db,p \ 

Sb.p 

c'^P-Xb.p 
\ c'^^y Yb^p j 



Cb' 



(2.11) 



where A'' is the number of unit cells, and the two phases in the right-hand side of the 
equation read = \{Px — t^) and ipy = \{py — tt). This transformation will be used 
in the next subsection for deriving the interaction Hamiltonian. 



2.2. The Heitler- London & Schubin-Wonsowsky-Zener interaction 

The Heitler-London (HL) interaction Hamiltonian describes the (intra- and inter- 
atomic) two-electron exchange. It comprises four parts [22, 23] corresponding to 
Cu 4s ^ O 2pcr, O 2pcr ^ Cu 3d^2_y2, O 2p^ ^ O 2py, and Cu 3da;2_j,2 ^ Cu 4s 
exchanges with transition amplitudes Jgp, Jp^, Jpp, and Jgd, respectively: 

n.ap 



'Jpd ^ ^Dl^aXl,pDnf3Xna + Dl^cXl^l3^nl3Yna 



n,af3 



-J, 



t A>t 



'^'^{n^,ny + l}a^nf:)-^{n:,.ny+l)l3Yna 
~^-^lna,,ny + l)a^('n^+l,ny)f3-^{n^,ny + l)0^{n^+l,ny}a 



(2.12) 



Let us now analyze the structure of the total electron Hamiltonian Ti.' = Ti-'^n + ''^hl ■ 
In terms of the Fermi operators 5"^, corresponding to the atomic orbitals, W reads: 



n 
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- E JvMp^J^c^jP^ (2-13) 

where iji = tjie^'^^', tji — t.^ and (f>ji = — (pi is the phase difference between 
the zth and jth atomic orbitals in the overlapping domain. Roughly speaking, onto 
every single-electron hopping amplitude one can map a corresponding two-electron 
hopping amplitude Jy. The case of a strong electron correlation implies that could 
be of the order of tij. Thus, one can expect that the following inequalities hold true 

In fact, the s-d exchange is the basic process responsible for the magnetism of 
transition metals; see for example reference [20]. It was understood since the dawn 
of quantum physics that the mechanism of ferromagnetism [30] is the two-electron 
exchange owing to the electron correlations [31]. 

Here we shall add a few words in retrospect concerning the two-electron correlation 
parameterized by in (2.13). Probably the first two-electron problem was Bohr's 
consideration of the He atom [32] (cf references [33, 34]) in which two electrons 
have opposite coordinates r2 = — Ti and momenta P2 = — Pi- For a purely radial 
motion, such a fall to the nucleus is stable and many years after Bohr's prediction 
double Rydberg states, with an effective Ry^g = (2 — 1/4) Ry, were discovered by 
electron energy loss spectroscopy [35]. These double Rydberg states with opposite 
electron momenta can be considered as proto-forms of the Cooper pairs. Interestingly, 
in 1914, Sir J J Thomson proposed [36] (cf. also the textbook [37]) that electric 
charge can propagate as electron doublets — another proto-form of the local (Ogg- 
Schafroth) pairs [38,39]. Before the appearance of quantum mechanics, Lewis [40] and 
Langmuir [41] introduced the idea of electron doublets in order to explain the nature 
of the chemical bond. Nearly at the same time Parson [42] came to the conclusion that 
"an electron is not merely an electron charge but a small magnet" or in his terminology 
"a magneton", cf. reference [40]. Later, in 1926, Lewis introduced also the notion of a 
photon [43] without any reliable theoretical background at the time. 

In the era of new quantum mechanics Heitler and London [21] realized the idea of 
electron doublets [44] and convincingly demonstrated how the two-particle correlation 
owing to a strong Coulomb repulsion can lead to a decrease of the energy, and by virtue 
of the Hellmann-Feynman theorem, to inter-atomic attraction for the singlet state of 
the electron doublet. The original Hcitler-London calculation, which is nowadays 
interpreted in every textbook in quantum mechanics and/or quantum chemistry, gives 
indeed a wrong sign of the exchange energy for very large inter-atomic distances but, 
in principle, there are no conceptual difficulties in the Heitler-London theory. The 
exchange energy J was represented [45] as a surface integral in the two-electron six- 
dimensional space (ri, r2) and this was shown to be an asymptotically exact result, cf. 
also reference [46]. The surface integral method gives amazingly accurate results (cf. 
the excellent monograph by Patil and Tang [47]) even if the exchange energy is of the 
order of the energies typical for solid state phenomena. Unfortunately, this method, 
that ought to be applied to ab initio calculated (e.g., from density functional theory 
(DFT) [48]) wave functions, is barely known in the solid state community (although 
a very recent work by Gor'kov and Krotkov [49] indicates that it is not completely 
forgotten) . 

This is one of the reasons why the t and J transfer integrals have been treated 
phenomenologically just as fitting parameters of the theory. A valuable discussion on 
a similar scope of ideas has recently been given by Brovetto, Maxia and Sails [50] but 
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it may well not be the only case. In order to ease comparison of the HL Hamiltonian 
with the other types discussed in the search of a theory of high-Tc superconductivity 
we shall rewrite it in terms of spin variables. 

The grounds for our theory have been set first by Schubin and Wonsowsky 
and later in more clear notions and notation by Zener [20]. The s-d two-electron 
exchange is the intra-atomic version of the HL interaction. Both of those 4-fermion 
interactions due to Heitler-London & Schubin- Wonsowsky-Zener can in principle 
mediate superconductivity and magnetism. 

2.2.1. Spin variables Let us introduce the spin operator Si and particle number 
operator fii for each atomic orbital, 

s^=^^ = ^ ^0 ^ = (*1t, *a) , (2.14) 

where ctq = l2x2 and cr are the Pauli sigma matrices, and the first two formulae imply 
summation over the spin indices. Introducing also the spin exchange operator Pij , 

Pi!^a.^^Jl3 = i!^|}i>Jo^, 4 = E ^^/^ ^-^a *j73 , (2.15) 

a/3 

we can rewrite the HL Hamiltonian per bond as [51-53] 

-JT. ^p^..^^c.^JP = = 2 J ■ S, -f \n,n\ . (2.16) 

We should stress that in the t-J model the term cx hifij enters with negative 
sign [54,56]. Let us also provide the "mixed" representation: 

+ ky [-i^x^jl + + {nn - %i) 

= + k-^^^n + kz (".T - ".i) , (2-17) 

where = ^^^^'ji, and 5,;^+ = ^It'^ii = ^ -■ l^ote that (2.16) imphes a purely 
orbital motion without spin flip: two electrons exchange their orbitals and only the 
spin indices reflect this correlated hopping. For J > 0, the HL Hamiltonian has a 
singlet ground state 

\S) = ^{i>\^i>]^-i>\^i>]^)\Y&c), 4',„|vac)=0, (vac|vac) = l, (2.18) 

with eigenvalue — J. The lowering in energy of the singlet state, having a symmetric 
orbital wave function, is of purely kinetic origin related to the delocalization of the 
particles at different orbitals. Symbolically, the "location" of the (approximately) 
localized electron doublet (s) in the structure signature of a molecule is designated by 
a colon, e.g.. H:H for the H2 molecule. This Lewis notation for the valence bond 
with energy — J (or four-Fermion terms in the second quantization language) is an 
important ingredient of the chemical intuition. In principle, such an exchange lowering 
is expected to exist for Bose particles as well. For electrons, however, we have triplet 
excited states 

|T+l) = *|^*]^|vac), 

|T0) =l^{i,\^i>]^ + i,m^)\y^c), (2.19) 
|T-1) =^^^^\ys.c), 
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with eigenvalue J. In the present work we consider the parameter J to be 
positive if it corresponds to antiferromagnetism, or pairing in the singlet channel. 
Thus the singlet-triplet splitting for the single-bond HL Hamiltonian (2.16) is 
2J. Similarly, the bonding-antibonding splitting for the single-particle hopping 
Hamiltonian — ^ X)a(^lQ^'ct + ^L^i") ^t, and the energy threshold for creation 
of a pair of normal carriers, considered in the next section, is 2A. Besides stemming 
from bare inter- and intra-atomic processes, two-electron hopping amplitudes J can 
be created by strong correlations [55] within the Hubbard model. For a nice review 
on this subject the reader is referred to the work by Spalek and Honig [56]. 



3. Reduced Hamiltonian and the BCS gap equation 

Substituting the Fermi operators ^'na, equation (2.11), into equation (2.12) one 
obtains the HL interaction Hamiltonian in a diagonal band representation. For the 
case of zero electric current [57], solely the reduced Hamiltonian TYhl-r, including 
creation and annihilation operators with opposite momenta only, has to be taken into 
account: 

^HL-R = ^ X] X] X] 4,p/34 -pa ^b,p;fc',p' h'-p'aCb',p'l3- (3.1) 
b.p b'.p' ap 

For singlet superconductors it is necessary to take into account the pairing with 
opposite spins, thereby the total reduced Hamiltonian reads 

''^R = X Vb,pcl,pa^b,pa + X X ^b,p;b',p' C^pi^ -piC6',-p'iC6',p'T , (3-2) 
b,p,a b.p b'.p' 

where rj^ p = tb,p ^ 1^ are the band energies measured from the chemical potential [57] . 
Hence the BCS equation [3] for the superconducting gap takes the familiar form 

^b,p = 77 y^(-Vb.p;b',p') ^'^ Afc/^p>, (3.3) 

jV f--' ^J^b'.p' 
b'.p' 

where E^^p ~ {Vh p + |A(),pl^)^^^ are the quasiparticle energies and n^.p = 
[exp(ii'b^p/fcgr) 4- 1]^"'^ are the Fermi filling factors with the Boltzmann constant, 
and T the temperature. The summation over the band index 6' should be restricted to 
the partially filled (metallic) bands, comprising sheets of the Fermi surface. Applying 
this standard procedure to the HL Hamiltonian (2.12), and after some algebra, we 
obtain the desired BCS pairing kernel, 

^b.p:b'.p' 2(/s(J SpSp' DpDp' '^pp ^x-^p-^p' ^y^p^p' 

+ 2 [JspSpSp' + JpdDpDp,) {-/xXpXp, + jyYpYp, ) , (3.4) 



. , fPy+P'y\ /q 

7^ = 4cos , 7y = 4cos — tt^ ] ■ (3-5) 



where 



As the band indices b and b' enter implicitly in the band energies e^^p in the equation 
for the eigenvectors ^p{eb.p); we will suppress them hereafter. The layered cuprates, 
admittedly, have a single conduction band and their Fermi surface has the shape of 
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a rounded square. In this simplest case one has to solve numerically the nonlinear 
integral equation 

— TT —TT 

The solution to this general gap equation, depending on the values, can exhibit s, 
p, or d-type symmetry. It has been shown previously that a purely p-p model [23, 58] 
( Jpp > 0) results in a <ixy {B2g) gap anisotropy. However, we found that an agreement 
with the experimentally observed dx2-y2 (Big) gap anisotropy (for a review see for 
example reference [59] ) can be achieved only in the simplest possible case of a dominant 
s-d exchange. This separable Hamiltonian deserves special attention and we will 
analyze it in the next sections. 



4. Separable s-d model 



For the special case of a purely s-d model, Jgp — Jpd = Jpp — 0, representing the 
spin exchange operator P as a (4 x 4)-matrix, cf. equations (2.15)-(2.16), the reduced 
pairing Hamiltonian takes the form 



p,q 







t 


/ 1 










s. 


-Pi Dpi 










1 





s. 


-p[Dpl 







1 










-plDpi / 




I 








1 / 



V S-qlbql 



(4.1) 



Carrying out an additional reduction for a spin-singlet pairing, the interaction 
Hamiltonian reads 



N 



(4.2) 



where —a stands for the electron spin projection opposite to a. For comparison, we 
provide again the kinetic energy part of the Hamiltonian employing the same notation. 



BH 



E 



/ Cd - M 





^pd^x 


ipd^y ^ 









- /" 


^sp^x 


^sp^y 






^pd-^a: 




ep - M 


^pp^x^y 




-^p,a 


\ ^tpdSy 


^sp^y 


^pp-^a; Sy 


Ep - M / 




V ^p,a / 



(4.3) 



Within the s-d model considered, the pairing kernel (3.4) factors into functions 
depending only on p or q. 



2 Jsd SpDp SqDq = 2 Jsd Xp Xq- 



(4.4) 



A schematic representation of the Jgd exchange amplitude is given in figure 2. 
This factorizable Markowitz-Kadanoff [60] form of the pairing kernel is a direct 
consequence of the local intra-atomic character of the s-d exchange in the transition 
ion. Substituting in equation (3.6) 

Ap{T) = S(T) SpDp = S(T) Xp, (4.5) 

one obtains in a closed form, cf. reference [61], a simple BCS equation for the 
temperature dependence of the gap, 

' xl , f E, ' 



2J, 



sd 



2En 



tanh 



p 

2KT , 



= 1, 



(4.6) 
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12 3 

r(A) 

Figure 2. Pairing two-clcctron exchange amplitude Jg^j "hidden" in the Cu atom, 
(a) Classical Bohr-Sommerfeld representation of the s-d two-electron exchange 
process. The inset shows how the Coulomb scattering leads to an effective electron 
spin exchange, (b) Electron charge distribution for Cu 4s and Cu 3d orbitals: the 
dashed line marks the Cu-O distance in the Cu02 plane. 



where 

s(,J + Aj)"'=[(«,-£p)^ + (E(T),Y,n"', (") 

Ep = /i. We wish to mention that separabihty of the order parameter (4.5) has 
been derived by Pokrovsky [62] in the general weak-couphng case and not only for 
factorizable pairing kernels. 

According to (2.4) we have 

Xp = SpDp = Asp tsp tpd {x - y) [es £p - 4ep tl^ {x + y) + 32ipp rf^ xy] 
X I [4ep tsp tpd {x -y)f + [e, 4 - 4ep tl^ {x + y) + 32<pp r.^^ xy] ^ 

+Ax [{Es £p - Sr^p y) ipd] ^ + 4y [(e, Sp - Sr^^ x) tpd] ^ } ^ (4.9) 

The gap symmetry is then easily made obvious in the narrow-band approximation. 
Formally, it is the asymptotic behaviour of the eigenvector (2.4) for vanishing hopping 
integrals t ^ Q. In this limit case [29], we have e3,p Cd, and 





1 


1 


\ 






"(^spipd/SsEp) ^ Sy 


) 


ATs.p 




(tpd/Sp) Sx 






V 


(ipd/Sp) Sy 


/ 



(4.10) 
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Figure 3. Electronic properties of the superconducting Cu02 plane. (a) 
Conduction band energy Ep as a function of the quasi-momentum p. The red 
contour corresponding to the Fermi energy, ep = Ep, is in excellent agreement 
with the ARPES data [63]. (b) Quasiparticle velocity function of 

quasimomcntum. The velocity variation along the Fermi contour is less than 10 %. 
The energy parameters are fitted to be in agreement with the typical ab initio 
calculations [64] . The significant overestimate disappears if the bandwidth is fitted 
to the experimental data, but the shape is conserved, (c) Momentum dependence 
of the gap-anisotropy function Xp within the s-d model. The functional values 
along the Fermi contour are indicated by a green line, (d) Superconducting gap at 
zero temperature Ap (green line) according to our analytical result (4.5), plotted 
along the Fermi contour (red line). The ARPES data [65] for Bi2Sr2CaCu208-j_i 
are given as prisms with sizes corresponding to the experimental error bars. The 
gap function along the Fermi contour has the same qualitative behavior and 
symmetry as the Cu 3d^2_j^2 electron wave function along the circular orbit 
sketched in figure 2(a). 



Clearly, D^ p exhibits Aig symmetry, while S^^p has the Big symmetry. Whence the 
product Ss^pD^ p oc cosp^, — cospy "inherits" the Big symmetry, figure 3(b), which is 
conserved even for realistic values of the hopping integrals, and from (4.5) it follows 

Ap cx Ss.pDs.p w — ?h2^ r(cosp^ - cospy). (4.11) 

(£/F — Csj (,-C/F — Cpj 

As can be seen in figure 3(c) this small-t approximation fits the ARPES data for the 
gap anisotropy quite well. Similar experimental data have been previously reported, 
e.g.; in reference [66]. Note, additionally, that close to the (tt, 7r)-point, where 
(Px — tt)^ + {py — tt)^ <C 1, the angular dependence of the gap can be written in 
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(4.12) 



The d-type angular dependence of both the gap anisotropy and the separable 
pairing kernel is often postulated in phenomenological model Hamiltonians to describe 
high-Tc superconductivity. The previous discussion thus provides a microscopic 
justification based on the fundamental exchange amplitudes. For the oxygen 
scenario [29,67], in which the Fermi level falls in a nonbonding oxygen band, 62. p ~ ep 
and t ^ [22], the gap function has different or additional nodes along the Fermi 
contour, 



S2,p 



4) 



(4.13) 



Here I?2.p and 6*2. p exhibit and symmetries, respectively. Let us also mention 
that the s- vector components, equation (2.3), constitute the arguments of the basis 
functions of the symmetry representations. 

Employing the analytical expression (2.8) for the constant-energy contours (CEC), 
one can implement an efficient numerical integration, 

2tt 2tt ct 



where 



dpx Jdpy f{ep) = Jde /(e), 

[{Ay + B^xjl -x) + (Ax + g)2y(l - y)] 



dp 



\A'xy + B'{x + y)+C'\ 



with A' , B' and C being the energy derivatives of the polynomials (2.9), 



-4'(e)-16[2t2ptpp_2t2^tpp 
+2 ^ 1.2 



ipp(£d + < 



= -4(<?p£d + <^d£.) - 4ep(i^p + <^d), 



(4.14) 



(4.15) 



(4.16) 



C'(e) =e,£2 _^ed£p + 2edesep• 
Using these functions, the band spectrum, see (2.8), can be obtained by Newton 



iterations 



Axy + B{x + y) + C 
A'xy^B'{x + y)+C' 



(4.17) 



with initial approximation for the conduction band £3 p = Ed- 

The charge carrier velocity is Vpao/fi, clq is the lattice constant, pi the 
dimensionless momentum component along the CEC, and Ch and et are the bottom 
and the top of the conduction band, respectively, eb < < ft- The canonic equation 
for the CEC (2.8), 



^ cos Pa, COSpj^ — {A+ 2B){cOSpj: 

can be cast in an explicit form 



cospy)+A + AC = 0, 



(4.18) 



PyAPx) = 2arcsin 



Bx + C 
Ax + B' 



PvAPx) = 27r - Pyjipx), (4.19) 
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and for the length element dpi we obtain 



The contour integration along the hole pocket Sp ~ const centered at the (tt, tt) point 
needs to be performed only over one eight of the CEC 

^dpi/(p.,p,) = 8 y'"/(p.,p,(p.)) (^^^) dp., (4.21) 



where 



Xd = sin^ , Axl + 2Bxd + C = 0. (4.22) 



5. Antiferromagnetic character of Jgd 

Let us address now the atomic physics underlying the s-d pairing mechanism. Within 
the framework of the Hartree-Fock (HF) theory the exchange energy is given [68] as 
an integral of the Cu 4s and Cu 3da;-2_y2 atomic wave functions, 

- -^id"^' = // rsiriWM ^^^ Mri)A{r2) dridr^, (5.1) 

and its sign corresponds to repulsion and depairing for singlet Cooper pairs. Thus, one 
can formulate the following conceptual problem, emerging in fundamental physics: 

(i) is it possible, as in the case of a covalcnt bond, for two-electron correlations to 
trigger a change of the sign of the exchange amplitude? 

(ii) how can one adapt the Heitler-London idea to a transition ion perturbed by 
ligands? 

There is no doubt that the solution to this problem will illuminate other problems in 
the physics of magnetism as well. In brief, the enigma can be stated as to whether 
the Heitler-London approximation for the exchange energy may result in Jgd > 0, cf. 
reference [68]. Let us recall that already in 1962 Herring [45] was advocating that 
"antiferromagnetic J^'s should be the rule, ferromagnetic J^-'s the exception". For 
the present, we can adopt the s-d model as a convenient microscopic phenomenology 
of superconductivity in the Cu02 plane. On the other hand, the exchange amplitude 
Jsd is an important ingredient in the physics of magnetism as well. 

Physics of magnetism certainly displays lots of subtleties, but for a qualitative 
comparison let us trace the "operation" of the s-d exchange amplitude Jsd in the case of 
the simplest model for a ferromagnetic metal. While for the Cu02 plane the s-band is 
empty, for transition metals it is the widest conduction band. The width of the d-band 
is significantly smaller and thus, making a caricature of the ferromagnetic metals, we 
completely neglect the width of the d-band. In this "heavy Fermion" approximation 
the d-electrons are considered as localized, and without significant energy loss they can 
be completely spin polarized, (ndt) ~ 1, {S'd,z) ~ 5 > 0- this case the self-consistent 
approximation applied to (2.17) gives 

2Sd ■ Ss ~ Sd,z {fisi - fisi) ~ (Sd^z) ("-sT - "si) ■ (5.2) 

Here Uga = (fisa) denotes the average number of s-electrons per atom with spin 
projection a. In order to calculate these variables one has to take into account the 



Superconductivity of overdoped cuprates 



14 



different filling of the s-bands with different polarizations, and sum over the quasi- 
momcnta. Finally, the exchange energy per atom reads 

Ex = -ijsd(»^sj, - »^st) < 0. (5.3) 

In the Cu02 plane, positive values of the Jgd parameter lead to singlet 
superconductivity. For ferromagnetic metals, positive values of Jgd correspond to 
polarization of the s-band opposite to d-state polarization, rigx — risT > 0. Thereby 
ferromagnetism could be brought about by an exchange amplitude with a sign 
corresponding to antiparallel spin polarization of s- and d-orbitals, cf. figure 4-15 
of reference [52]. Thus the same sign of the s-d exchange amplitude Jgd can be at 
the origin of ferromagnetism, e.g., in Fe and Ni, and superconductivity in the Cu02 
plane. This is perhaps the simplest scenario for cuprate superconductivity based on 
the two-electron exchange processes. 

According to a naive interpretation of Hund's rule the Kondo effect should not 
exist. In the epoch-making paper [69] on the resistance minimum in dilute magnetic 
alloys Kondo concluded that in the s-d exchange model, due to Zener [20] , Kasuya and 
Yosida [70], the sign of the direct exchange amplitude Jgd must be antiferromagnetic. 
And vice versa, the minimum disappears if Jgd is ferromagnetic. Such minimum 
exists for many magnetic metals and alloys and is another hint in favor of Herring's 
argument [45] mentioned earlier. In his analysis Kondo speculates that Jgd is a 
parameter whose sign and magnitude have to be determined so as to fit the experiment, 
and concluded that antiferromagnetic values of the order of eV are quite reasonable. 
For a review on the Kondo problem we refer the reader to reference [71]. 

On the other hand every textbook on atomic physics tells us that parallel electron 
spins and an antisymmetric wave function minimize the electrostatic energy. Put 
differently, the tendency toward ferromagnetism in Hund's rule is of electrostatic 
origin. As Kondo has pointed out [69], the problem is to find the origin of an 
antiferromagnetic Jgd, or how to overcome the strong electrostatic repulsion. It is 
very plausible that it is not a single driving force, but instead one has to take into 
account several interfering electron scattering amplitudes. 

5.1. Intra-atomic correlations 

The self-consistent approximation has been known in celestial mechanics for ages. 
Accordingly, the motion of a planet is averaged over its orbit. One has to calculate 
then the potential created by this orbital-averaged motion and perform a sum over all 
particles. Where does this scheme fail? It fails in the case of a resonance when the 
periods for some planets are commensurate or just equal. This is nothing but the case 
of a transition ion for which the energies and classical periods are very close. Then 
the resonant repetitive electron scattering, symbolically presented in figure 2, leads to 
strong electron correlations like in the double Rydberg states of atoms [35] . For double 
Rydberg states in He it is necessary to solve a two-electron quantum problem but for 
other atoms we have to take into account the influence of the other electrons in some 
self-consistent approximation, the local density approximation (LDA), for example. 
For two 4s electrons the two-electron correlations are so strong that they have to 
be taken into account from the very beginning [35]. There are no doubts that the 
two-electron correlations between 4s and 3d electrons having almost equal energies 
cannot be neglected. In other words Hartree-Fock theory cannot be used directly. 
Hence, the Bohr picture is not merely a historical remark but rather an indispensable 
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ingredient of tlic contemporary physics of magnetism: two-electron correlations can be 
important even in a single atom. We thus conclude that the two-electron correlations 
may overrule Hund's rule for the local s-d exchange. Note also that the single-particle 
orbitals (accessible, e.g., from DFT [48], HF and Xa methods, etc.) only form an 
adequate basis for a subsequent account of electronic correlations. A first step in this 
direction will be ab initio calculation of the Jed- 

5.2. Indirect s-d exchange 

The antiferromagnetism of the insulating phase of the undoped cuprates is mediated by 
the Bloch-Kramers- Anderson indirect exchange [72] between 3d electrons of nearest- 
neighbour Cu ions via O 2p electrons. It is unlikely that the numerical value of this 
Jdd exchange integral be dramatically changed in the metallic phase obtained by hole 
doping. In the metallic phase, however, the same indirect exchange mechanism will 
operate between 3d and 4s electrons at the same Cu atom via the 2p electrons of the 
O ligands. For illustration, let us compare the indirect s-d exchange amplitude t/g^"'^' 
with Jdd- There are three important factors: (i) Every Cu ion has four O ligands, 
figure 1 (a) . (ii) The hopping amplitude between 4s and 2p orbitals is bigger than the 
3d-2p transfer, (iii) The Cu on-site Coulomb repulsion between 4s and 3d electrons 
Usd is much smaller than the 3d-3d Hubbard repulsion Udd- Taking into account these 
factors one can expect that J^™'^^ is an order of magnitude bigger than Jdd- 



The relatively small Jdd ensures Neel temperatures Tn of the order of room 
temperature. Hence we can conclude that the indirect exchange can contribute 
significantly to the total Jgd amplitude responsible for the pairing. However, only 
very detailed first-principles calculations can clarify the relative contributions of the 
direct and the indirect s-d exchange. 

5.3. Effect of mixing wave functions 

In an early paper [73], by analyzing the 5-shift and the anomalous Hall effect in Gd 
metal, Kondo showed that an antiferromagnetic Jgd can result from the effect of mixing 
the wave functions of conduction and d electrons. We believe that this property is 
preserved if the d electrons also form a conduction band, or even in the case of a single 
s-p-d hybridized band. We should note that Kondo's argumentation for the need of 
a Jsd with an antiferromagnetic sign in the Kondo effect is related to Anderson's 
consideration of localized magnetic states in metals [74]. In the latter schematized 
model, based on the works of Friedel [75] , Anderson shows that "any g-shift caused by 
free-electron polarization will tend to have antiferromagnetic sign." 

As it was expected by various investigators the later numerical calculations 
confirmed that the striking features of negative hyperfine field with large amplitude 
comes mainly from the contact contribution of the core electrons [76-78]. The 
antipolarization between the s- and d-electrons in transition metal compounds is also 
well observed by Mossbauer spectroscopy, however the contribution of the core s- 
electrons and conduction band cannot be experimentally resolved. For the pairing, 
the amplitude of the s-d Kondo scattering is essential because in some terminology 
the Cu02 plane can be considered to be a single-band Kondo lattice, cf. reference [79]. 




(5.4) 
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Given the above diversity of channels for s-d exchange it is not surprising that 
an adequate first-principles scheme to calculate Jgd is still sought. Furthermore, Jgd 
is involved in the theory of magnetism in an entangled way precluding so far a direct 
relation between ab initio calculations and formulae fitting the experiment [80]. 

5.4- Cooper and Kondo singlet formation 

Although Kondo [69] does cite Zener's paper [20], in many publications in the field of 
magnetism the s-d interaction is referred to as Kondo interaction. Often this term is 
used in a broader sense causing eventually terminological misunderstanding. Here we 
cite some works which could be related to our theory even though the relation may not 
be direct or immediately apparent. Analyzing the possible "interplay of Cooper and 
Kondo singlet formations in high- Tc cuprates" Sekitani et al [81] point out that "In 
the 20^^ century, two significant many-body phenomena due to spin singlet formation 
were discovered in the field of solid state physics: superconductivity and the Kondo 
effect". These authors believe that the pseudogap in the normal state corresponds to 
the dissociation energy of the Kondo bound states and that superconductivity and the 
Kondo effect are competing in Lai.85Sro.i5Cu04. They speculate that the interplay 
between superconductivity and the Kondo effect has not be considered seriously for 
high-Tc cuprates (further references on the Kondo effect in cuprate superconductors 
can be found, e.g., in [82]). It would be premature for us to comment on these ideas; 
we note however that within this terminology our theory could be considered as a 
Kondo interaction mechanism for pairing in high-Tc cuprates. 

It is unclear from microscopic point of view if the same "Kondo interaction" 
amplitude Jgd is responsible for the empirical Kadowaki- Woods ratio [83] but the 
location of Lai.7Sro.3Cu04 on the A (the coefficient of the resistivity) versus 70 
(the electronic specific-hcat coefficient) plot, cf. figure 4 in [84], is a serious hint that 
Jsd in the Cu02 plane is one of the largest exchange amplitudes in solid state physics, 
comparable with the uranium heavy- fermion compounds and SrVOs. In this sense 
our theory requires a large, yet acceptable Jsd value, putting the cuprates among the 
most interesting materials with considerable exchange interaction. 

6. Dogmatics, Discussion, Conclusions and Perspectives 

In a review on the history of studies of superconductivity and the prospects for further 
research in the field Ginzburg [2] conditionally divided the history into several periods: 

(i) The "Day Before Yesterday" (1911-1941). This period starts with the discovery of 
superconductivity in Leiden by Gilles Hoist and Heike Kamerlingh Onnes [2,85]. 

(ii) "Yesterday" (1942-1986). This period embraces the appearance of the '^-/S. 
theories and the first significant technical applications. 

(iii) "Today" (1987-?). This epoch emerged with the discovery of the high-Tc 
cuprates [1]. 

(iv) "Tomorrow" (?). The final landmark of "Today" must be some event. 

Long ago, in the "Day Before Yesterday" high-Tc superconductivity was known as 
a "foZue dream" of physicists. Considerable theoretical efforts were applied "Yesterday", 
attempting to predict possible realizations of this phenomenon [86]. At that time the 
problem of high-Tc superconductivity was "one of the most interesting and attractive 
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problems from the purely scientific point of view" [86] . Intriguingly, the special role 
of layered metallic systems and almost two-dimensional superconductivity [87] was 
mentioned already in this epoch, and a big variety of mechanisms of superconductivity 
were considered including the s-d exchange [88,89]. This exchange process was well 
known in the physics of magnetism since the dawn of quantum mechanics. Thus it 
is not surprising that the first work on the s-d pairing mechanism, by Akhiezer and 
Pomeranchuk [88], was accomplished about a year after the celebrated BCS paper [3]. 
These pioneering works, however, "have been ignored thus far" [90]. 

After Bednorz and Miiller's work [1] the problem of high-Tc superconductivity 
soon came into fashion. After experiencing the 'smell of roast meat yesterday skeptics 
or even critics can become zealous advocates of a new direction of endeavor. But this 
is another story — more in the realm of psychology and sociology than scientific and 
technical activity" [91]. All models of high-Tc superconductivity were revisited in 
great detail in the uncountable number of papers that have appeared in the epoch 
"Today". 

6.1. Aesthetics and frustrations of the central dogmas 

The common trends of some new theoretical models for cuprate superconductivity 
were systematized by Anderson [92] in six dogmas. We find it very instructive to 
compare our theory of high-Tc superconductivity with these dogmas. 

"Dogma I: All the relevant carriers of both spin and electricity reside in the Cu02 
plane and derive from the hybridized 2p-Cu Sd^^-y^ orbital which dominates the 
binding in these compounds. ...in summary look at the planes only (a great and 
welcome simplification.) " 

The key ingredient of our pairing theory is the four-fermion s-d interaction 
between the Cu 4s and Cu 3dj,2_j,2 orbitals. If we cut the Cu 4s orbital off from 
the Hilbert space of the Cu02 plane such a pairing interaction cannot exist. Although 
Cu 4s is an empty band, it is an important component of the theory of high-Tc 
superconductivity. The O 2p orbitals are the intermediaries between the Cu 4s and 
Cu 3d2.2_j,2 orbitals that create the necessary s-d hybridization of the conduction 
Cu 'id^2_y2 band. 

"Dogma II: Magnetism and high-Tc superconductivity are closely related, in a very 
specific sense: i.e., the electrons which exhibit magnetism are the same as the charge 
carriers. . . . We must solve the old problem of doping of a single Mott-Hubbard band 
before we can begin the problem of high-Tc." 

The incommensurate spin-density waves (SDW) observed in the superconducting 
phase of La2CuO4.11 and La1.88Sro.12O4 by neutron scattering [93] and muon spin 
relaxation [94], respectively, demonstrated that antiferromagnetism of the Cu site is 
innocuous for superconductivity in the cuprates. These antiferromagnetic correlations 
are not depairing and do not change significantly Tc and the electronic structure 
of Cu02 plane. The observed correspondence between the magnetic and the 
superconducting order parameters is an additional hint that both phenomena have 
a common origin [95]; see also the detailed theoretical works [96]. Nevertheless 
the coexistence of SDW and superconductivity with a common critical temperature 
cannot be clearly observed in every high-Tc cuprate. As a result superconductivity 
can be considered, at least in first approximation, separately from a possible 
antiferromagnetism as is done in the present paper. In Cr metal the amplitude of 
the SDW shows also a BCS-like temperature dependence [97] and the SDW-theory is 
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based on the conventional theory of metals based in turn on the Landau Fermi-liquid 
theory. We consider the quasiparticlc picture as a reliable starting point for the theory 
of high-Tc cuprates as well. 

"Dogma III: The dominant interactions are repulsive and their energy scales are 
all large. . . . Restrict your attention to a single band, repulsive (not too big) U Hubbard 
Model. " 

Indeed, the dominant interactions are repulsive — "Nobody has abrogated the 
Coulomb law", as Landau used to emphasize [2]. However, something subtle occurs 
when the atomic orbitals are analyzed. The strong electron repulsion leads to Heitler- 
London type correlations: two electrons cannot occupy simultaneously the same 
orbital, not even if they have opposite spins. The exchange of electrons between 
two orbitals decreases the electron kinetic energy and thereby the total energy of the 
whole system. In molecular physics, according to the Hellmann-Feynman theorem such 
a decrease in energy drives an inter-atomic attraction for large inter-atomic distances. 
Thus, the valence attraction is the final result of the dominant Coulomb repulsion 
between electrons. In this way the Heitler-London-type exchange between itinerant 
electrons gives rise to electron-electron attraction and conventional Cooper pairing. 
The s-d exchange, "residing" in the Cu atom, can be considered as an "intra-atomic- 
valence bond" — an attraction-sign scattering amplitude due to the Coulomb repulsion 
between the correlated electrons. The s-d exchange in the transition ions is one of the 
most intensive exchange processes in solid state physics. Such a high-frequency process 
is described by the exchange amplitude Jsd in the lattice models for the electronic 
structure and its sign is determined by the inter-electronic Coulomb repulsion. The 
Heitler-London interaction is a result of strong electron repulsion and survives even for 
infinite Hubbard U. This interaction is lost when starting with the infinite-C/ Hubbard 
model, however. Thus, not a single-band Hubbard model but a single-band s-d model 
with antiferromagnetic exchange amplitude is the adequate starting point for a realistic 
treatment of Cu02 superconductivity. 

"Dogma IV: The 'normal' metal above T^. . . is not a Fermi liquid. . . but retains a 
Fermi surface satisfying Luttinger's theorem at least in the highest-Tc materials. We 
call this a Luttinger Liquid. " 

Very recently, the crucial experiment has finally been conducted. After 15 years 
of intensive investigations of the cuprates it is now experimentally established [99] 
that the overdoped cuprates obey the 150-years-old Wiedemann-Franz law within 
a remarkable 1% accuracy. After this experimental clarification the theoretical 
comprehension will hardly keep us waiting long. This experiment has also solved 
the old problem of the nature of charge carriers created by doping of a single Mott- 
Hubbard band, cf. Dogma II. Now we know that charge carriers of the normal state are 
standard Landau quasi-particles [100] for which we have conventional Cooper pairing 
in the superconducting phase. "Holons", "spinous" and spin-charge separation are 
unlikely to occur and behave so as to emulate the properties of the ideal Fermi gas. 
As a function of the hole doping per Cu atom, p, the critical temperature is a smooth 
parabola [98], 

Tjj-nii,^ w 1 - 82.6(p - 0.16)^ (6.1) 

Thus, it is improbable that the nature of the carriers and pairing mechanism can be 
dramatically changed in the optimal and underdoped regimes although a number of 
new and interesting phenomena complicate the physics of the underdoped cuprates. 
In short, in our opinion the experimental validation of the Wiedemann-Franz 
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law in overdoped cuprates [99] is a triumph of the Landau [100] and Migdal concept 
of Fermi quasi-particles (and Landau spirit of trivialism in general) and provides a 
refutation of the spin-charge separation in cuprates [101]. Hence, the problem of 
deriving the Wiedemann-Franz law for strongly correlated electrons in the Cu02 
plane has just been set in the agenda. According to the Fermi liquid theory [102] 
interactions between the particles create an effective self-consistent Hamiltonian. As 
Kadanoff [103] has pointed out, this idea was much developed by Landau [104] and 
Anderson [105]. Unfortunately, for high-Tc cuprates a link is still missing between the 
Landau quasiparticle concept and the one due to Slater that even scattering matrix 
elements can be calculated from first principles. 

"Dogma V: Nonetheless, enough directions have been probed to indicate strongly 
that this odd-even splitting of Cu02 planar states doesn't exist. . . . The impact of 
Dogma V, then, is that the two-dimensional state has separation of charge and spin 
into excitations which are meaningful only within their two-dimensional substrate; to 
hop coherently as an electron to another plane is not possible, since the electron is a 
composite object, not an elementary excitation. " 

Within the single-particle approximation (section 2) the bilayer band splitting is 
readily obtained from (2.8) and (4.19) by the replacements 

ei^e,±t±^ii, i^s,p,d, (6.2) 

where t_L,ii is the hopping amplitude between the ith orbitals in the adjacent Cu02 
planes. In other words, the two constant energy curves due to the bilayer splitting are 
described by the same equation (2.8). Since it is plausible that t±^ss dominates, from 
(4.10) one finds 

Ai^biiayor ~ St^^^s I "Ss.p | ^ ~ 22 mcV (cosp^ - cospy)^ (6.3) 

in agreement with references [64,106]. The numerical value of 22 meV has been 
reported for heavily overdoped Bi2Sr2CaCu208+5 (BSCCO) [106]. This experiment, 
crucial for Dogma V, cf. reference [107], is another piece of evidence in favor of the 
conventional behaviour of the electron excitations in the (Cu02)2 slab. Since Ai?biiaycr 
is relatively small in comparison with the width of the conduction band, it is another 
hint that even for bilayer superconductors like BSCCO and YBa2Cu307_5 (YBCO) 
the analysis of a single Cu02 plane is an acceptable initial approximation. Similarly, 
for fitting the three-dimensional Fermi surface of Tl2Ba2Cu06+5 determined by angle 
magnetoresistance oscillations [108] one can start with the simplest possible tight- 
binding approximation 

£i £( + t_i__ii cospz, i ^ s, p, d. (6.4) 

"Dogma VI: Interlay er hopping together with the "confinement" of Dogma V 
is either the mechanism of or at least a major contributor to superconducting 
condensation energy. " 

The interlayer hopping which is understood as a single-electron process definitely 
cannot be considered as a two-electron pairing interaction creating the condensation 
energy. It is only one of the details when one concentrates on the material-specific 
effects in high-Tc superconductors. The inter-slab hopping between double (Cu02)2 
layers is a coherent Josephson tunnelling responsible for the long-living plasma 
oscillations with frequency Wpi < A. These plasma oscillations along with far infrared 
transparency of the superconducting phase were theoretically predicted [109] for 
BSCCO — one of the few predictions made for high-Tc cuprates, cf. the postdiction 
[110]. After the experimental observation [111], the plasma resonances associated 
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with the Cooper-pair motion soon turned into a broad research field [112]. Subgap 
plasmons were predicted [113] for conventional superconducting thin films as well, and 
shortly after experimentally confirmed [114] for thin Al films on SrTiOa substrate. 
The relatively lagged development of the physics of this effect was partially due to the 
false neglect of the longitudinal current response in the classical works on microscopic 
theory. Concluding, let us note that the London penetration depth A can be considered 
as the Compton wave length of the Higgs boson of mass mnc^ = Swpi, but the overall 
contribution of the interlaycr hopping to the condensation energy is negligible. 

6.2. Discussion 

The band structure of the Cu02 plane is now believed to be understood. However, 
after 15 years of development a mismatch of a factor of two or three between 
the ab initio and the experimental spectroscopic estimates for the single-electron 
hopping amplitudes t, or the bandwidth, tends to be interpreted rather as a state- 
of-the-art "coincidence". The Hcitler-London approach is well-known in quantum 
chemistry [115, 116], and has been successfully used for a long time in the physics 
of magnetism [117]. We hope that realistic first-principles calculations aiming at the 
exchange integrals J of the Cu02 plane can be easily carried out. Should they validate 
the correct (antifcrromagnetic) sign and the correct order of magnitude of Jsd, we can 
consider the theory of high-Tc superconductivity established. We stress that the two- 
electron exchange, analyzed here, is completely different from the double exchange 
considered in reference [118]. 

In order to compare the derived results with the experiment, it is necessary that 
the tight-binding conduction band energy be fitted to the available ARPES data. In 
doing so a few parameters have to be properly taken into account: the Fermi energy 
Ep, as determined from the total area of the hole Fermi contour, the difference between 
the Fermi energy and the Van Hove singularity, Ep — e(7r,0), the difference between 
the Van Hove singular point and the bottom of the conduction band at the F point, 
e(7r,0) — e(0,0). The fit may further allow taking into account a possible realization 
of the Abrikosov-Falkovsky scenario, cf. reference [29]. According to the latter, for 
Ed < Ep < Ss and sufficiently small tpd, the conduction band can be the narrow 
(nonbonding) oxygen band having a perfect (within the framework of the four-band 
model) extended Van Hove singularity. If the superconducting gap has a _Big-type 
symmetry, its maximum value along the Fermi contour. A^ax — max]Ap(T — 0)], 
determines the Jsd exchange integral in the s-d model. Thus, the temperature 
dependence of the gap, described by the hmction S(r), and the overall thermodynamic 
behaviour and low frequency electrodynamic response will be determined without free 
fitting parameters. 

The fit to the extended Van Hove singularity as observed, e.g., in reference [119] 
also points to a relatively low-lying Cu 4s level and one needs to consider the minimum 
value of Eg- Although the 4s band is completely empty (Cu 4s level is above the Fermi 
level), a very close location is not "harmless" and would necessarily lead to some 
prediction for the optical behaviour. With some risk of opening the Pandora's box, 
we should mention that the lowest position of the Cu 4s level is determined by the 
mid-infrared (MIR) response. According to this possible interpretation, the broadly 
discussed maximum of the absorption in the MIR range is due to 3d-4s interband 
transition: one electron in the conduction band is excited by the light to the empty 
Cu 4s band. It seems that, up to now, there is no natural explanation of this MIR 
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optical adsorbtion (for a review see [8]). 

The derived gap anisotropy function (4.9) and its interpolation (4.11) compared to 
the ARPES experiment showed that the "standard" four-band model spanned on the 
Cu 3dx2_y2, Cu 4s, O 2p2;, and O 2py orbitals, with an antifcrromagnctic s-d pairing 
interaction, successfully describes the main features of the ARPES data: the rounded- 
square-shaped Fermi surface, small energy dispersion along the (0, 0)-(27r, 0) line, and 
the d-typc (Big) symmetry of the energy gap Ap along the Fermi contour. According 
to the pairing scenario proposed here, strong electron correlations "drive" the electron 
exchange amplitudes. These inter- and intra-atomic processes occur on energy 
scales unusually large for solid state physics. However, the subsequent treatment 
of the lattice Hamiltonian can be performed completely within the framework of the 
traditional BCS theory. The criterion for applicability of the BCS scheme is not given 
by the J vs t, but rather by the Tc vs Ep — eh relation. Taking into account the 
typical ARPES-derived bandwidths, which are much bigger than Tc we come to the 
conclusion that the BCS trial wave function [120] is applicable for the description of 
superconductivity in the layered cuprates with an acceptable accuracy if Tc does not 
significantly exceed room temperature. 

It is worth adding also a few remarks on the normal properties of the layered 
cuprates. Among all debated issues in the complex physics of the cuprates, the 
most important one is perhaps that of the normal-phase kinetics. The long-standing 
problem is whether the paring interaction dominates (or totally determines) the 
mechanism of Ohmic resistance in the normal phase, as is the case for conventional 
superconductors. Within the present theory this question can be formulated as follows: 
does the s-d exchange interaction dominate the scattering of the normal-state charge 
carriers above Tc? This is a solvable kinetic problem whose rigorous treatment will 
be given elsewhere. Here we shall restrain ourselves in providing only a qualitative 
discussion. 

In electron-electron scattering, just like in traffic accidents, the crucial effect comes 
from the backscattering in "head-on" collisions. For backscattering (i.e., = tt) in the 
case of s-d interaction, it turns out that the matrix elements entering the pairing 
amplitude are also important. It can be easily realized that this amplitude vanishes 
along the diagonals of the Brillouin zone (0, 0)-(7r, tt) (the F-M direction). Thereby the 
cold spots on the Fermi contour correspond to the zeros of Ap. And vice versa, the hot 
spots are associated with a maximum gap along (tt, 7r)~(0, tt) (the M-X direction). In 
this sense, cuprates repeat the qualitative feature of the conventional superconductors, 
with a maximal gap corresponding to maximal scattering on the Fermi surface. 

All layered cuprates are strongly anisotropic and two-dimensional models give a 
reasonable starting point to analyze the related electronic processes. Most importantly, 
the picture of a layered metal brings in something qualitatively new which does 
not exist for a bulk metal — the "interstitial" electric field between the layers, like 
the one in any plane capacitor. The thermodynamic fluctuations of this electric 
field and related fluctuations of the electric potential and charge density constitute 
an intensive scattering mechanism analogous to the blue-sky mechanism of light 
scattering by density fluctuations. It has recently been demonstrated [121] that the 
experimentally observed linear resistance can be rationalized in terms of the plane 
capacitor scenario; density fluctuations in the layered conductors are more important 
than the nature of the interaction. In such a way the linear normal-state resistivity 
is an intrinsic property [121] of the "layered" electron gas and cannot be used as an 
argument in favor of non-Fermi-liquid behavior. The resistance of the normal phase 
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may not be directly related to the pairing mechanism and these problems can be 
solved separately. Nevertheless it will be interesting to check whether the anisotropic 
scattering in cuprates [122-124] can be explained within the framework of the s-d 
pairing Hamiltonian. 

The present theory can also predict a significant isotope effect in the cuprates. 
Even though the Jsd pairing amplitude does not depend on the atomic mass, the 
charge carriers reside the ionic Cu02 2D lattice, thereby rendering polaron effects, as 
in any ionic crystal, possible. For the lighter oxygen isotope the lattice polarization 
is more pronounced, leading to enhanced effective mass and density of states, and 
reducing the transfer integrals. Overall, the isotope effect in the Cu02 plane is due 
to the isotope effect of the density of states. This rationale can be quantitatively 
substantiated. In our theory, upon isotope substitution at T = 0, e.g., the change of 
the penetration depth would be mainly driven by the Cooper pair effective mass that 
could be determined by means of the Bernoulli effect. At that the superfluid density 
remains unchanged. The calculation of the isotope effect on Tc requires an evaluation 
of the polaron effects on the conduction band. Although this is a feasible problem, it 
is beyond the scope of the present work. 

The proposed mechanism for pairing in the Cu02 plane can be handled 
much like an "Alice-in- Wonderland" toy-model, but we find it fascinating that all 
ingredients of our theory are achievements of quantum mechanics dating back to the 
memorable 1920s, presently described in every physical textbook, and constituting 
the fundamentals of solid state physics [125, 126]. It would be worthwhile attempting 
to apply the approach, used in this paper, for modelling triplet and heavy-fermion 
superconductivity as well. 

6.3. Conclusions: the reason for the success of the Cu02 plane 

We find it very instructive to analyze qualitatively the reasons for the success of the 
realization of high-Tc superconductivity in the Cu02 plane: 

(i) Because of the relatively narrow quasi-two-dimensional conduction d-band, due 
to p-d hybridization, the density of states is rather high. The wide s-band 
resulting from s-p hybridization is completely empty, which is somewhat unusual 
for compounds containing transition ions. 

(ii) The pairing s-d exchange process was known since the first years of quantum 
physics. It is omnipresent in the physics of the transition ions but in order for it 
to become the pairing mechanism in perovskites it is necessary that the s- and 
d-levels be close. In other words, a virtual population of the s-level is at least 
needed in order to make the Jsd amplitude operative. Indeed, the conduction d- 
band is, actually, a result of the s-p-d hybridization in the two-dimensional Cu02 
plane. 

With the above remarks, one can speculate that among the perovskites the layered 
ones are more favorable for achieving higher Tc(cf. also the discussion in the Appendix 
on page 26). The transition ion scries ends with Cu'^^ and the Cu 3(1,^2 _y2 and 4s levels 
are too close. One should keep in mind that the filling of the electron shells finishes 
with a "robbery" in Cu [52]: 3d-'^°4s^ instead of 3d^4s^ as one could expect from the 
electron configuration of the Ni atom (3d^4s^). However, the energy difference between 
these two Cu shell configurations is very small. 
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Another favorable factor is the proximity of the O 2p and Cu 3d^2_j,2 levels. 
Thus, post factum the success of Cu and O looks quite deterministic: the Cu02 
plane is a tool to realize a narrow d-band with a strong s-p-d hybridization. It was 
mentioned earlier that Jsd is one of the largest exchange amplitudes, but the 4s and 
3d orbitals are orthogonal and necessarily require an intermediary whose role is played 
by the O 2p orbital. Hence this theory can be nicknamed "the 3d-to-4s-by-2p highway 
to superconductivity^^ [19]. The Jsd amplitude is omnipresent for all transition ion 
compounds, the hybridization of 3d, 4s and 2p is however specific only for the Cu02 
plane. 

How this qualitative picture can be employed to predict new superconducting 
compounds is difficult to assess immediately. We believe, however, that this picture, 
working well for the overdoped regime, is robust enough against the inclusion of all 
the accessories inherent to the physics of optimally doped and underdoped cuprates: 
cohabitation of superconductivity and magnetism [95], stripes [127], pseudo-gap [128], 
interplay of magnetism and superconductivity at individual impurity atoms [129], apex 
oxygen, Cu02 plane dimpling, doping in chains [130], the 41 meV resonance [131], 
etc. Perhaps some of these ingredients can be used in the analysis of triplet 
superconductivity in the copper-free layered perovskite Sr2Ru04 [132]. It is also likely 
that the superconductivity of the Ru02 plane is a manifestation of a ferromagnetic 
exchange integral J. The two-electron exchange mediates superconductivity and 
magnetism in heavy Fermion compounds [133] as well. We suppose that lattice models 
similar to the approach here will be of use in revealing the electronic processes in these 
interesting materials. Two-electron exchange may even contribute to the 30 K Tc of 
the cubic perovskite Bao.eKo. 46103 but so far it is difficult to separate the exchange 
contribution from the phonon part of the pairing interaction. However, the strange 
doping behaviour of Tl2Ba2Cu06±5 in comparison with YBCO requires more detailed 
investigation [11]. 

6.4- Perspectives: if "Tomorrow" comes. . . 

The technological success in preparing the second generation of high-Tc superconduct- 
ing cables by depositing thin-layer superconducting ceramics on a flexible low-cost 
metallic substrate is crucial for the future energy applications. The USA Department 
of Energy suggests global superconducting energy products would command an annual 
market of 30 G$ by about 2020. High-Tc superconductor power cables, transformers, 
motors and generators could grab a 50% market share by 2013, 2015, 2016 and 2021, 
respectively [134]. On the other hand atomic-layer engineering of superconducting 
oxides will trigger progress in materials science and electronics. One can envision 
multi-functional all-oxide electronics, e.g., sensors, processing and memory devices, 
all monolithically integrated within a single chip [135]. In spite of the technological 
progress and tens of thousands of publications the theoretical "picture in early 2000 
remains fairly cloudy on the whole'''' [2]. The landmark of "Today" must be some 
event. "What event will it be? It is desirable that this landmark be the insight into 
the mechanism of superconductivity in high-Tc cuprates" [2] . 

In this paper we presented a traditional theory for superconductivity in overdoped, 
and possibly also optimally doped cuprates. All of its ingredients can be found 
in the textbooks and there is a considerable chance that we witness the victory 
of traditionalism, as it was in the history of quantum electrodynamics (QED) 
half a century ago, but it may well be just a personal viewpoint "brainwashed by 
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Feynman" [136]. Nonetheless let us use the example of QED to illustrate the essence 
of our contribution. QED appeared as a synthesis between perturbation theory 
and relativity. Both components had been known well before the QED conception. 
Similarly, both the BCS theory and the exchange interaction have been known for 
ages, so the point in the agenda was how to conceive out of them the theory 
of high-Tc cuprates. Such a theory contains necessarily a big number of energy 
parameters {Ef, Cs, ep, Cd, tsp, tpp, tpd, Jsd, Jpd, Jsp, Jpp) which are difficult to 
determine simultaneously! (for the current status of the problem see for example 
reference [137,138]). The first step will definitely be to use ARPES data in which 
the spectrum is clearly seen and to neglect in a first approximation the "irrelevant" 
inter-atomic exchange integrals Jpd, Jsp and Jpp. In this case, for a known normal 
spectrum one can determine Jsd from Tc or from the maximum gap at T = 0. 

A crucial "meeting point" between theory and experiment is the Ginzburg- 
Landau (GL) theory. The general form of the GL coefficients for anisotropic-gap 
superconductors, including the effect of disorder, is given in reference [139] and is 
directly applicable to the present model. For the s-d separable kernel (4.4) the specific 
heat C{T) in the clean limit can also be explicitly derived [140] and has the GL form 
C{T) = Cn + Ca, with 

^?N(r) =y(g,(z.p)). 



Ca{T) ^k^T— ^ ———^ — -, — --e{Tc~T), 6.5) 
b 7C(3 Xp96 I'p) 



where Vp = jk^, and 



^ ^(^MM), h{T) = _m_(4,,(,^)), (6.6) 
, . 1 , , TT^ 1 / tanh V 1 \ Q v"^ 

^'^y^i = — = tttt^^ tit- ' 9c(i^) 



2 cosh^ V ' 14C(3) v'^ \ V cosh^ v) ' tt^ cosh^ v ' 

Accordingly, the jump of the specific heat at Tc is expressed by the GL coefficients a 
and 6, AC = k^T^a^{^^)lh{^^. With the help of the general equations (6.6) one can 
further determine the influence of the Van Hove singularity on the thermodynamic 
and electrodynamic behaviour. For AC the effect of the Van Hove singularity is 
reported in reference [141] , and for a general review on the Van Hove scenario of high- 
Tc superconductivity we refer the reader to reference [142]. When the Fermi level is 
not close to the Van Hove singularity the GL coefficients can be worked out as integrals 
over the Fermi surface; methodological details are given in reference [143]. Knowledge 
of the GL coefficients is also fundamental for the physics of fluctuation phenomena in 
superconductors [144]. 

Furthermore, a microscopic consideration of the London penetration depth A for 
screening currents in the Cu02 plane gives [145] 

1 / dp; 

~'ord{vp)——^, (6.7) 



A2(r) eoC^^i^rfeff/ ^ ^'(27r)2' 

X The gap-anisotropy fit in figure 3(d) is quite robust against the choice of the parameters. To 
illustrate and emphasize the capability of the model we have used, for example, unrealistically big 
values of the hopping integrals: = 0, £s = 5, £p = —0.9, tpd = 1.13, tsp = 1.63, tpp = 0.2 eV. This 
set of parameters corresponds to band calculations but gives a factor 2—3 wider conduction band. If 
the band is fitted to the ARPES data Jg^ can be less than 1 eV. A realistic fit is deemed to be a 
subject of a collaboration with experimentalists. 
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ri=0 

where the integration is performed along the Fermi contour. The penetration depth 
A(T) is involved in the Bernoulli effect in superconductors [146]: 



7^H 2eoc2 ^ ' 7^H ag^ff 7^''^"'(2^)2 
where Aiy9 is the change of the electric potential induced by a current density j, 
1/7?.H = eritot is the volume charge density of the charge carriers, with doff the effective 
spacing between the Cu02 planes. 

For given penetration depth extrapolated to zero temperature, A(0), and Hall 
constant of the superconducting phase, one can easily determine the effective mass of 
the Cooper pairs 

^* = -^^ |e1=2|e|. (6.9) 

This important material parameter m* is experimentally accessible from the 
electrostatic modulation of the kinetic inductance of thin superconducting films [147] 
as well as from the surface Hall effect [148]. 

Having a big variety of calculated variables the parameters of the theory can 
be reliably fitted. Another research direction is the first-principles calculation of 
the transfer amplitudes and two-electron exchange integrals. The level of agreement 
with the fitted values will be indicative for the completeness of our understanding. 
In addressing more realistic problems, the properties of a single space-homogeneous 
Cu02 plane will be a reasonable starting point. Concluding, we believe that there is a 
true perspective for the theoretical physics of cuprate superconductors to become an 
important ingredient of their materials science. 

Magnetism and superconductivity are among the most important collective 
phenomena in condensed matter physics. And, remarkably, magnetism of transition 
metals and high-Tc superconductivity of cuprates seem to be two faces of the same 
ubiquitous two-electron exchange amplitude. 
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Appendix. Tc-eg correlations: a hint toward the mechanism of 
superconductivity in cuprates 

In the letter by Pavarini et al , reference [149], a strong correlation is observed between 
Tcmax and a single parameter 

s(e) = (6,-6)(6-ep)/(2t,p)2, 

which is controlled by the energy of the Cu 4s orbital eg, Cp is the O 2p single cite 
energy, ts-p is the transfer integral between neighbor Cu 4s and O 2p orbitals, and 
Cp is the conduction band energy whose bottom at the F point coincides with the Cu 
Zd.x2_y2 energy ed- It is unfortunate that theorists have not so far paid any attention to 
this observation because it is an important correlation between the ah initio calculated 
parameter r ~ (1 + s)/2 and the experimentally measured Tc max which can reveal the 
subtle link between the experiment and the theory and finally solve the long-standing 
puzzle of the mechanism of high-Tc superconductivity (HTSC). 

The purpose of the present Comment [150] is to emphasize that the missing link 
has already been found, and the work by Pavarini et al can be used as a crucial 
test for theoretical models of HTSC. Perhaps the simplest possible interpretation, 
though one could search for alternatives, is given within the framework of the present 
theory. In order for the Schubin-Zener-Kondo exchange amplitude Jsd to operate as 
a pairing interaction of the charge carriers, it is necessary the Cu 4s orbital to be 
significantly hybridized with the conduction band. The degree of this hybridization 
depends strongly on the proximity of the Cu 4s level to the Fermi level ey. Thus, it 
is not surprising that eg controls the maximal critical temperature max; being the 
only parameter of the Cu02 plane which is essentially changed for different cuprate 
superconductors . 

Cu 3d and Cu 4s are orthogonal orbitals and their hybridization is indirect. 
First, the Cu 3da;2_j,2 orbital hybridizes with the O 2pa; and O 2pj^ orbitals, then 
the O 2p orbitals hybridize with Cu 4s. As a result we have a "3d-to-4s-by-2p" 
hybridization of the conduction band of HTSC cuprates which makes it possible 
the strong antiferromagnetic amplitude Jsd to create pairing in a relatively narrow 
Cu 3da;2_j,2 conduction band. The hybridization "filling" of the Cu 4s orbitals can 
be seen in cluster calculations as well [151]. The s(e) parameter introduced in 
reference [149] reflects the proximity of all 3 levels in the generic 4-band Hamiltonian 
of the Cu02 plane — their "random coincidence" for the Cu-0 combination. Suppose 
that those levels are not so close to each other. In this case the slightly modified 
parameter 

s'(6)-(eg-e)(e-ep)/(4tgpipd) (.1) 

is simply the energy denominator of the perturbation theory which describes the 
hybridization filling of the axial orbital, see (4.10). Whence [s'(e)]^ is a denominator 
of the pairing amplitude in the BCS equation (6.4). Hence, we conclude that the 
correlations reported in reference [149] are simply the correlations between the critical 
temperature Tc and the dimensionless BCS coupling constant p(eF) Jsd/[s'(eF)]^. Of 
course, for coupling constants ~ 1 the BCS trial wavefunction can be used only for 
qualitative estimates, but knowing the Hamiltonian the mathematical problem may 
somehow be solved. In any case, even qualitatively, we are sure that the stronger 
pairing amplitude Jsd and hybridization l/s'^ep) enhance Tc- 

Having LDA calculations for the band structures of many cuprates it is worthwhile 
performing a LCAO fit to them [64] and using experimental values of to extract 
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the pairing amplitude Jgd for all those compounds. The ab initio calculation of the 
Kondo scattering amplitude parameterized by Jgd is an important problem which has 
to be set in the agenda of computational solid state physics. We expect that it will 
be a weakly material dependent parameter of the order of the s-d exchange amplitude 
in Kondo alloys, but perhaps slightly bigger as for the Cu ion the 3d and 4s levels 
are closer compared to many other ions. Closer energy levels, from classical point of 
view, imply closer classical periods of orbital motion which leads to some "resonance" 
enhancement of the exchange amplitude due to intra-atomic two-electron correlations. 
The final qualitative conclusion that can be extracted from the correlations reported 
by Pavarini et al is the explanation why only the Cu02 plane renders HTSC possible, 
whereas hundreds other similar compounds are not even superconducting, or have 
only a "conventional" value of Tc- The natural explanation is: because its s-parameter 
is not small enough below its critical value. Even among the cuprates one can find 
compounds with "conventional" values of Tc having relatively large value of the s- 
parameter. For other transition metal compounds the parameter 

s'(ed) = (cs - ed)(ed - ep)/(4tspipd) 

is much bigger than its critical value Sc which can be reached probably only for the 
Cu-0 combination. Thereby, the correlation reported by Pavarini et al is a crucial 
hint which of the models for HTSC is still on the arena. Concluding, we also note that 
the 4s hybridization is responsible for the three-dimensional coherent Fermi surface of 
TlaBazCuOe+a [152] 
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